Isotopic signatures induced by upwelling reveal regional fish stocks in Lake Tanganyika.

Lake Tanganyika's pelagic fish sustain the second largest inland fishery in Africa and are under pressure from heavy fishing and global warming related increases in stratification. The strength of water column stratification varies regionally, with a more stratified north and an upwelling-driven, biologically more productive south. Only little is known about whether such regional hydrodynamic regimes induce ecological or genetic differences among populations of highly mobile, pelagic fish inhabiting these different areas. Here, we examine whether the regional contrasts leave distinct isotopic imprints in the pelagic fish of Lake Tanganyika, which may reveal differences in diet or lipid content. We conducted two lake-wide campaigns during different seasons and collected physical, nutrient, chlorophyll, phytoplankton and zooplankton data. Additionally, we analyzed the pelagic fish-the clupeids Stolothrissa tanganicae, Limnothrissa miodon and four Lates species-for their isotopic and elemental carbon (C) and nitrogen (N) compositions. The δ13C values were significantly higher in the productive south after the upwelling/mixing period across all trophic levels, implying that the fish have regional foraging grounds, and thus record these latitudinal isotope gradients. By combining our isotope data with previous genetic results showing little geographic structure, we demonstrate that the fish reside in a region for a season or longer. Between specimens from the north and south we found no strong evidence for varying trophic levels or lipid contents, based on their bulk δ15N and C:N ratios. We suggest that the development of regional trophic or physiological differences may be inhibited by the lake-wide gene flow on the long term. Overall, our findings show that the pelagic fish species, despite not showing evidence for genetic structure at the basin scale, form regional stocks at the seasonal timescales. This implies that sustainable management strategies may consider adopting regional fishing quotas.

Reference genome choice and filtering thresholds jointly influence phylogenomic analyses.

Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant-calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate to what extent the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find greatest topological accuracy when filtering sites for minor allele count >3-4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with minor allele count >1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short read genomic data for phylogenetic inference.

Cycles of fusion and fission enabled rapid parallel adaptive radiations in African cichlids.

Although some lineages of animals and plants have made impressive adaptive radiations when provided with ecological opportunity, the propensities to radiate vary profoundly among lineages for unknown reasons. In Africa's Lake Victoria region, one cichlid lineage radiated in every lake, with the largest radiation taking place in a lake less than 16,000 years old. We show that all of its ecological guilds evolved in situ. Cycles of lineage fusion through admixture and lineage fission through speciation characterize the history of the radiation. It was jump-started when several swamp-dwelling refugial populations, each of which were of older hybrid descent, met in the newly forming lake, where they fused into a single population, resuspending old admixture variation. Each population contributed a different set of ancient alleles from which a new adaptive radiation assembled in record time, involving additional fusion-fission cycles. We argue that repeated fusion-fission cycles in the history of a lineage make adaptive radiation fast and predictable.

Conceptual and empirical bridges between micro- and macroevolution.

Explaining broad molecular, phenotypic and species biodiversity patterns necessitates a unifying framework spanning multiple evolutionary scales. Here we argue that although substantial effort has been made to reconcile microevolution and macroevolution, much work remains to identify the links between biological processes at play. We highlight four major questions of evolutionary biology whose solutions require conceptual bridges between micro and macroevolution. We review potential avenues for future research to establish how mechanisms at one scale (drift, mutation, migration, selection) translate to processes at the other scale (speciation, extinction, biogeographic dispersal) and vice versa. We propose ways in which current comparative methods to infer molecular evolution, phenotypic evolution and species diversification could be improved to specifically address these questions. We conclude that researchers are in a better position than ever before to build a synthesis to understand how microevolutionary dynamics unfold over millions of years.

Parallel shifts in trout feeding morphology suggest rapid adaptation to alpine lake environments.

Eco-evolutionary interactions following ecosystem change provide critical insight into the ability of organisms to adapt to shifting resource landscapes. Here we explore evidence for the rapid parallel evolution of trout feeding morphology following eco-evolutionary interactions with zooplankton in alpine lakes stocked at different points in time in the Wind River Range (Wyoming, USA). In this system, trout predation has altered the zooplankton species community and driven a decrease in average zooplankton size. In some lakes that were stocked decades ago, we find shifts in gill raker traits consistent with the hypothesis that trout have rapidly adapted to exploit available smaller-bodied zooplankton more effectively. We explore this morphological response in multiple lake populations across two species of trout (cutthroat trout, Oncorhynchus clarkii, and golden trout Oncorhynchus aguabonita) and examine the impact of resource availability on morphological variation in gill raker number among lakes. Furthermore, we present genetic data to provide evidence that historically stocked cutthroat trout populations likely derive from multiple population sources, and incorporate variation from genomic relatedness in our exploration of environmental predictors of feeding morphology. These findings describe rapid adaptation and eco-evolutionary interactions in trout and document an evolutionary response to novel, contemporary ecosystem change.

Complete mitochondrial genomes and updated divergence time of the two freshwater clupeids endemic to Lake Tanganyika (Africa) suggest intralacustrine speciation.

BACKGROUND: The hydrogeological history of Lake Tanganyika paints a complex image of several colonization and adaptive radiation events. The initial basin was formed around 9-12 million years ago (MYA) from the predecessor of the Malagarasi-Congo River and only 5-6 MYA, its sub-basins fused to produce the clear, deep waters of today. Next to the well-known radiations of cichlid fishes, the lake also harbours a modest clade of only two clupeid species, Stolothrissa tanganicae and Limnothrissa miodon. They are members of Pellonulini, a tribe of clupeid fishes that mostly occur in freshwater and that colonized West and Central-Africa during a period of high sea levels during the Cenozoic. There is no consensus on the phylogenetic relationships between members of Pellonulini and the timing of the colonization of Lake Tanganyika by clupeids. RESULTS: We use short-read next generation sequencing of 10X Chromium libraries to sequence and assemble the full mitochondrial genomes of S. tanganicae and L. miodon. We then use Maximum likelihood and Bayesian inference to place them into the phylogeny of Pellonulini and other clupeiforms, taking advantage of all available full mitochondrial clupeiform genomes. We identify Potamothrissa obtusirostris as the closest living relative of the Tanganyika sardines and confirm paraphyly for Microthrissa. We estimate the divergence of the Tanganyika sardines around 3.64 MYA [95% CI: 0.99, 6.29], and from P. obtusirostris around 10.92 MYA [95% CI: 6.37-15.48]. CONCLUSIONS: These estimates imply that the ancestor of the Tanganyika sardines diverged from a riverine ancestor and entered the proto-lake Tanganyika around the time of its formation from the Malagarasi-Congo River, and diverged into the two extant species at the onset of deep clearwater conditions. Our results prompt a more thorough examination of the relationships within Pellonulini, and the new mitochondrial genomes provide an important resource for the future study of this tribe, e.g. as a reference for species identification, genetic diversity, and macroevolutionary studies.

Variable hybridization between two Lake Tanganyikan cichlid species in recent secondary contact.

Closely related taxa frequently exist in sympatry before the evolution of robust reproductive barriers, which can lead to substantial gene flow. Post-divergence gene flow can promote several disparate trajectories of divergence ranging from the erosion of distinctiveness and eventual collapse of the taxa to the strengthening of reproductive isolation. Among many relevant factors, understanding the demographic history of divergence (e.g. divergence time and extent of historical gene flow) can be particularly informative when examining contemporary gene flow between closely related taxa because this history can influence gene flow's prevalence and consequences. Here, we used genotyping-by-sequencing data to investigate speciation and contemporary hybridization in two closely related and sympatrically distributed Lake Tanganyikan cichlid species in the genus Petrochromis. Demographic modelling supported a speciation scenario involving divergence in isolation followed by secondary contact with bidirectional gene flow. Further investigation of this recent gene flow found evidence of ongoing hybridization between the species that varied in extent between different co-occurring populations. Relationships between abundance and the degree of admixture across populations suggest that the availability of conspecific mates may influence patterns of hybridization. These results, together with the observation that sets of recently diverged cichlid taxa are generally geographically separated in the lake, suggest that ongoing speciation in Lake Tanganyikan cichlids relies on initial spatial isolation. Additionally, the spatial heterogeneity of admixture between the Petrochromis species illustrates the complexities of hybridization when species are in recent secondary contact.

Hybridization decreases native cutthroat trout reproductive fitness.

Examining natural selection in wild populations is challenging, but crucial to understanding many ecological and evolutionary processes. Additionally, in hybridizing populations, natural selection may be an important determinant of the eventual outcome of hybridization. We characterized several components of relative fitness in hybridizing populations of Yellowstone cutthroat trout and rainbow trout in an effort to better understand the prolonged persistence of both parental species despite predictions of extirpation. Thousands of genomic loci enabled precise quantification of hybrid status in adult and subsequent juvenile generations; a subset of those data also identified parent-offspring relationships. We used linear models and simulations to assess the effects of ancestry on reproductive output and mate choice decisions. We found a relatively low number of late-stage (F3+) hybrids and an excess of F2 juveniles relative to the adult generation in one location, which suggests the presence of hybrid breakdown decreasing the fitness of F2+ hybrids later in life. Assessments of reproductive output showed that Yellowstone cutthroat trout are more likely to successfully reproduce and produce slightly more offspring than their rainbow trout and hybrid counterparts. Mate choice appeared to be largely random, though we did find statistical support for slight female preference for males of similar ancestry. Together, these results show that native Yellowstone cutthroat trout are able to outperform rainbow trout in terms of reproduction and suggest that management action to exclude rainbow trout from spawning locations may bolster the now-rare Yellowstone cutthroat trout.

The Genetic Population Structure of Lake Tanganyika's Lates Species Flock, an Endemic Radiation of Pelagic Top Predators.

Understanding genetic connectivity plays a crucial role in species conservation decisions, and genetic connectivity is an important component of modern fisheries management. In this study, we investigated the population genetics of four endemic Lates species of Lake Tanganyika (Lates stappersii, L. microlepis, L. mariae, and L. angustifrons) using reduced-representation genomic sequencing methods. We find the four species to be strongly differentiated from one another (mean interspecific FST = 0.665), with no evidence for contemporary admixture. We also find evidence for strong genetic structure within L. mariae, with the majority of individuals from the most southern sampling site forming a genetic group that is distinct from the individuals at other sampling sites. We find evidence for much weaker structure within the other three species (L. stappersii, L. microlepis, and L. angustifrons). Our ability to detect this weak structure despite small and unbalanced sample sizes and imprecise geographic sampling locations suggests the possibility for further structure undetected in our study. We call for further research into the origins of the genetic differentiation in these four species-particularly that of L. mariae-which may be important for conservation and management of this culturally and economically important clade of fishes.

Investigating the morphological and genetic divergence of arctic char (Salvelinus alpinus) populations in lakes of arctic Alaska.

Polymorphism facilitates coexistence of divergent morphs (e.g., phenotypes) of the same species by minimizing intraspecific competition, especially when resources are limiting. Arctic char (Salvelinus sp.) are a Holarctic fish often forming morphologically, and sometimes genetically, divergent morphs. In this study, we assessed the morphological and genetic diversity and divergence of 263 individuals from seven populations of arctic char with varying length-frequency distributions across two distinct groups of lakes in northern Alaska. Despite close geographic proximity, each lake group occurs on landscapes with different glacial ages and surface water connectivity, and thus was likely colonized by fishes at different times. Across lakes, a continuum of physical (e.g., lake area, maximum depth) and biological characteristics (e.g., primary productivity, fish density) exists, likely contributing to characteristics of present-day char populations. Although some lakes exhibit bimodal size distributions, using model-based clustering of morphometric traits corrected for allometry, we did not detect morphological differences within and across char populations. Genomic analyses using 15,934 SNPs obtained from genotyping by sequencing demonstrated differences among lake groups related to historical biogeography, but within lake groups and within individual lakes, genetic differentiation was not related to total body length. We used PERMANOVA to identify environmental and biological factors related to observed char size structure. Significant predictors included water transparency (i.e., a primary productivity proxy), char density (fish·ha-1), and lake group. Larger char occurred in lakes with greater primary production and lower char densities, suggesting less intraspecific competition and resource limitation. Thus, char populations in more productive and connected lakes may prove more stable to environmental changes, relative to food-limited and closed lakes, if lake productivity increases concomitantly. Our findings provide some of the first descriptions of genomic characteristics of char populations in arctic Alaska, and offer important consideration for the persistence of these populations for subsistence and conservation.

Structural genomic variation leads to genetic differentiation in Lake Tanganyika's sardines.

Identifying patterns in genetic structure and the genetic basis of ecological adaptation is a core goal of evolutionary biology and can inform the management and conservation of species that are vulnerable to population declines exacerbated by climate change. We used reduced-representation genomic sequencing methods to gain a better understanding of genetic structure among and within populations of Lake Tanganyika's two sardine species, Limnothrissa miodon and Stolothrissa tanganicae. Samples of these ecologically and economically important species were collected across the length of Lake Tanganyika, as well as from nearby Lake Kivu, where L. miodon was introduced in 1959. Our results reveal differentiation within both S. tanganicae and L. miodon that is not explained by geography. Instead, this genetic differentiation is due to the presence of large sex-specific regions in the genomes of both species, but involving different polymorphic sites in each species. Our results therefore indicate rapidly evolving XY sex determination in the two species. Additionally, we found evidence of a large chromosomal rearrangement in L. miodon, creating two homokaryotypes and one heterokaryotype. We found all karyotypes throughout Lake Tanganyika, but the frequencies vary along a north-south gradient and differ substantially in the introduced Lake Kivu population. We do not find evidence for significant isolation by distance, even over the hundreds of kilometres covered by our sampling, but we do find shallow population structure.

The coincidence of ecological opportunity with hybridization explains rapid adaptive radiation in Lake Mweru cichlid fishes.

The process of adaptive radiation was classically hypothesized to require isolation of a lineage from its source (no gene flow) and from related species (no competition). Alternatively, hybridization between species may generate genetic variation that facilitates adaptive radiation. Here we study haplochromine cichlid assemblages in two African Great Lakes to test these hypotheses. Greater biotic isolation (fewer lineages) predicts fewer constraints by competition and hence more ecological opportunity in Lake Bangweulu, whereas opportunity for hybridization predicts increased genetic potential in Lake Mweru. In Lake Bangweulu, we find no evidence for hybridization but also no adaptive radiation. We show that the Bangweulu lineages also colonized Lake Mweru, where they hybridized with Congolese lineages and then underwent multiple adaptive radiations that are strikingly complementary in ecology and morphology. Our data suggest that the presence of several related lineages does not necessarily prevent adaptive radiation, although it constrains the trajectories of morphological diversification. It might instead facilitate adaptive radiation when hybridization generates genetic variation, without which radiation may start much later, progress more slowly or never occur.

Playing out Liem's Paradox: Opportunistic Piscivory across Lake Tanganyikan Cichlids.

Trophic specialization is a key feature of the diversity of cichlid fish adaptive radiations. However, K. F. Liem observed that even species with highly specialized trophic morphologies have dietary flexibility, enabling them to exploit episodic food resources opportunistically. Evidence for dietary flexibility comes largely from laboratory studies, and it is unclear whether cichlid fishes undergo diet shifts in the wild. We report observations of diet switching by multiple cichlid species in Lake Tanganyika as a consequence of unusual concentrations of schooling juvenile clupeid fishes. Fish species with varying degrees of trophic specialization converged on a single prey: juvenile sardines that are also endemic to Lake Tanganyika (Stolothrissa tanganicae and Limnothrissa miodon). We provide evidence for cichlid species acting as jacks-of-all-trades and discuss this evidence in the framework of Liem's classic paradox: that trophic specialization does not preclude dietary flexibility.

Variable hybridization outcomes in trout are predicted by historical fish stocking and environmental context.

Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still do not know how variable outcomes of hybridization are across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co-occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while trout in other locations had a broader range of ancestry, including both parental species and first-generation hybrids. Later-generation intermediate hybrids were rare relative to backcrossed hybrids and rainbow trout individuals. Using an individual-based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under null expectations of random mating and no selection against hybrids. Since this deviation implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modelled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.

Rapid buildup of sympatric species diversity in Alpine whitefish.

Adaptive radiations in postglacial fish offer excellent settings to study the evolutionary mechanisms involved in the rapid buildup of sympatric species diversity from a single lineage. Here, we address this by exploring the genetic and ecological structure of the largest Alpine whitefish radiation known, that of Lakes Brienz and Thun, using microsatellite data of more than 2000 whitefish caught during extensive species-targeted and habitat-randomized fishing campaigns. We find six strongly genetically and ecologically differentiated species, four of which occur in both lakes, and one of which was previously unknown. These four exhibit clines of genetic differentiation that are paralleled in clines of eco-morphological and reproductive niche differentiation, consistent with models of sympatric ecological speciation along environmental gradients. In Lake Thun, we find two additional species, a profundal specialist and a species introduced in the 1930s from another Alpine whitefish radiation. Strong genetic differentiation between this introduced species and all native species of Lake Thun suggests that reproductive isolation can evolve among allopatric whitefish species within 15,000 years and persist in secondary sympatry. Consistent with speciation theory, we find stronger correlations between genetic and ecological differentiation for sympatrically than for allopatrically evolved species.

The smelly path to sympatric speciation?

In this issue of Molecular Ecology, Poelstra et al. () use genomic data to show that cichlid species in the adaptive radiation in Lake Ejagham, Cameroon, experienced gene flow with a riverine relative, primarily prior to their diversification. Intriguingly, this introgression brought with it olfactory alleles that the authors suggest may play a sensory role in speciation. As a classic example of sympatric speciation due to the highly restricted geography of this small (0.49 km2 ) crater lake (Schliewen et al., ), this result sheds new light on the history of this fascinating radiation. As genomic data for cichlid radiations accumulate, finding evidence of introgression increasingly appears to be the rule rather than the exception (Kautt et al., ; Meier et al., , ). However, two points are pressing as evidence for introgression mounts: (a) It is crucial to understand the assumptions of analytical approaches used and (b) it is important to think clearly about the divergent contexts in which evidence for introgression has been invoked as an important feature of adaptive radiation.

Divergent parasite infections in sympatric cichlid species in Lake Victoria.

Parasitism has been proposed as a factor in host speciation, as an agent affecting coexistence of host species in species-rich communities and as a driver of post-speciation diversification. Young adaptive radiations of closely related host species of varying ecological and genomic differentiation provide interesting opportunities to explore interactions between patterns of parasitism, divergence and coexistence of sympatric host species. Here, we explored patterns in ectoparasitism in a community of 16 fully sympatric cichlid species at Makobe Island in Lake Victoria, a model system of vertebrate adaptive radiation. We asked whether host niche, host abundance or host genetic differentiation explains variation in infection patterns. We found significant differences in infections, the magnitude of which was weakly correlated with the extent of genomic divergence between the host species, but more strongly with the main ecological gradient, water depth. These effects were most evident with infections of Cichlidogyrus monogeneans, whereas the only host species with a strictly crevice-dwelling niche, Pundamilia pundamilia, deviated from the general negative relationship between depth and parasitism. In accordance with the Janzen-Connell hypothesis, we also found that host abundance tended to be positively associated with infections in some parasite taxa. Data on the Pundamilia sister species pairs from three other islands with variable degrees of habitat (crevice) specialization suggested that the lower parasite abundance of P. pundamilia at Makobe could result from both habitat specialization and the evolution of specific resistance. Our results support influences of host genetic differentiation and host ecology in determining infections in this diverse community of sympatric cichlid species.

Pattern and Process in the Comparative Study of Convergent Evolution.

Understanding processes that have shaped broad-scale biodiversity patterns is a fundamental goal in evolutionary biology. The development of phylogenetic comparative methods has yielded a tool kit for analyzing contemporary patterns by explicitly modeling processes of change in the past, providing neontologists tools for asking questions previously accessible only for select taxa via the fossil record or laboratory experimentation. The comparative approach, however, differs operationally from alternative approaches to studying convergence in that, for studies of only extant species, convergence must be inferred using evolutionary process models rather than being directly measured. As a result, investigation of evolutionary pattern and process cannot be decoupled in comparative studies of convergence, even though such a decoupling could in theory guard against adaptationist bias. Assumptions about evolutionary process underlying comparative tools can shape the inference of convergent pattern in sometimes profound ways and can color interpretation of such patterns. We discuss these issues and other limitations common to most phylogenetic comparative approaches and suggest ways that they can be avoided in practice. We conclude by promoting a multipronged approach to studying convergence that integrates comparative methods with complementary tests of evolutionary mechanisms and includes ecological and biogeographical perspectives. Carefully employed, the comparative method remains a powerful tool for enriching our understanding of convergence in macroevolution, especially for investigation of why convergence occurs in some settings but not others.

Evolution in a Community Context: On Integrating Ecological Interactions and Macroevolution.

Despite a conceptual understanding that evolution and species interactions are inextricably linked, it remains challenging to study ecological and evolutionary dynamics together over long temporal scales. In this review, we argue that, despite inherent challenges associated with reconstructing historical processes, the interplay of ecology and evolution is central to our understanding of macroevolution and community coexistence, and cannot be safely ignored in community and comparative phylogenetic studies. We highlight new research avenues that foster greater consideration of both ecological and evolutionary dynamics as processes that occur along branches of phylogenetic trees. By promoting new ways forward using this perspective, we hope to inspire further integration that creatively co-utilizes phylogenies and ecological data to study eco-evolutionary dynamics over time and space.